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Abstract 

An effective simulation of quantum entanglement is presented using classical fields modulated 
with n pseudorandom phase sequences (PPSs) that constitute a n2”-dimensional Hilbert space with 
a tensor product structure. Applications to classical fields are examplied by effective simulation of 
both Bell and GHZ states, and a correlation analysis was performed to characterize the simulation. 
Results that strictly comply with criteria of quantum entanglement were obtained and the approach 
was also shown to be applicable to a system consisting of n quantum particles. 
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Quantum entanglement, one of the most fascinating and important features in quantum 
theory, is widely appreciated as an essential ingredient in quantum computations Ml. 
Simulations of quantum entanglement through optical approaches were investigated both 
theoretically and experimentally jb-l?]. A quantum bit can be represented by a distinct 
)ath or space mode of a classical held in an interferometric setup as classical optics analogies 
9l-ll7|. However, a n-qubit system with 2" basis states must be represented by 2” distinct 
paths or modes of a classical held. These simulations are usually not ehective due to an 
exponential increase in required physical resources correlated with the addition of quantum 
bits jo, 181. The drawback can be traced to a lack of a rigorous tensor-product structure of 


the system 


18 


19|. It is very inspiring that polarized beams (radially and azimuthally) 


of classical held exhibit a tensor product structure and are isomorphic to the Bell states by 
adding some degrees of freedom of a single system 201. 


In this letter, we present an ehective simulation of quantum entanglement of n quantum 
bits by using an analogy of classical helds modulated with pseudo-random phase sequences 
(BBSs). Based on the properties of BBSs, we proved that the n helds modulated with n 
diherent BBSs constitute a n2”-dimensional Hilbert space with a tensor product structure, 
which dihers signihcantly from those in classical simulations that were executed lacking a 
tensor-product structure jo i?]. By using an optical interferometric setup, BBSs yield not 
only random measurement results, but also an ensemble model to dehne the ensemble average 
and correlation functions 2l|. The BBSs, derived from orthogonal pseudorandom sequences. 


are widely applied to Code Division Mu 


a way to distinguish diherent users [22 


tiple Access (CDMA) communication technology as 


24| . A set of pseudorandom sequences is generated 


by using a shift register guided by a Galois held GF{p) that satishes orthogonal, closure 


and balance properties 


23| . In this letter, we utilize a m-sequence of period N — 1(A^ = 


generated by a primitive polynomial of degree s over GF{p) and apply it to binary phase 
shift modulation, a well-known modulation format in wireless and optical communications 
24|. Next we generate a BBS set S = over GF{2), where each A^*^ 

is a phase sequence with N phase units and time slots: A*^*^ = [Aj^*\ A^*^ • • • A^^^], while A*^^^ is 


22 


an all-0 sequence and other sequences can be generated by using following method 


23 


2 ^: 


(1) given a primitive polynomial of degree s over GF{2), a base sequence of a length 2^ — 1 is 
generated by the Linear Feedback Shift Register; (2) other sequences are obtained by cyclic 
shifting of the base sequence; (3) by adding zeroes to the sequences, the occurrence of any 


2 








































element equals to 2^ — 1; (4) mapping the elements of the sequences to [0,27r]: 0 mapping 
0, 1 mapping vr. 

We first consider two orthogonal modes (polarization or transverse modes), |0) and 
|1), of a classical held. A simulation state can be expressed as a mode superposition: 
IV’) = a |0) + /9|1), where jap + |/3p = 1, (a,/3 G C). All of the mode superposition 
states span a Hilbert space, where we will explore properties associated with this special 
classical held. By introducing a map / : A —)■ on the set of S, we obtain a phase sequence 
set H = l...A^}, and with which were written a superposed state 

corresponding to n-th sequence: 

|V-„) = e‘-"”’K|0)+/3„|l)) (1) 


According to the properties of m-sequence, the set H has following properties: (a) closure: 
the product of any two sequences equals one of sequence in the set; (b) balance: except 


any sequence of set H satisfy = 0,V6* G M; (c) orthogonality: 

any two of the sequences satisfy normalized correlation: 9^1*Vfc 

which equals 1 when i = j, and 0 otherwise. In fact, the map / corresponds to phase 
modulations of BBSs of S onto the classical held. 

A BBS map / constitutes a phase ensemble, wherein each phase unit represents a single 
simulation, and measurement of a physical quantity is a result of ensemble average. Similar 


to t 


lat in quantum mechanics, ensemble average and correlation measurement can be dehned 


22|-|24| . In the quadrature demodulation, each code obtained in a sequence unit of a BBS 


can be considered as a single measurement. The sequence number of the BBS’s unit can be 
used to label the ensemble. Diherent from the ergodicity hypothesis of quantum mechanics, 
the ergodicity of BBS is determined and much more efficient. 

Given the properties of the BBSs and the Hilbert space, the inner product of any two 
states and their orthogonal property can be obtained by: 


N 


k=l 



( 2 ) 


where A^*\ are the fc-th units of A^*) and respectively. Based on above properties, 
the classical helds modulated with different BBSs are independent and distinguishable. Fig. 
1 shows construction pathway of simulation states, generated by unitary transformed from 
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initial states-the mode superposition of classical fields with PPS Furthermore, a 

general form of a simulation state can be constructed from \'ipn)'- 


N 


N 


i=i 


(3) 


2=1 


Following the pathway in Fig. 1, a simulation state |\h) is obtained, denoting with a direct 
product of iV'n): 

I'h) = \'lpl) 0 \'lp 2 ) • • • 0 IV’iv) (4) 

Due to closure property of PPS, the phase sequence of each state 1*1*2 • • Dat) remains 
to belong to D and consists of the product of multiple sequences. Basis for Hilbert space of 
simulation is spanned by |*i *2 • • 0Af) |j = 1 • • • = Oorlj, with a total base state 

number of N2^. Generally a simulation state takes the form: 


21=0 2lV=0 


N 


Kp 2 • • Gjv) 

.i=i 


(5) 


where ,• denotes a total of N2^ coefficients. It is obvious that the Hilbert simulation 
space is greater than what is required for simulation of quantum state. To obtain a space 
the same size as that in quantum mechanics, either restrictions or a proper measurement 
need to apply j2l[ |. 

PPS provides not only the tensor structure and space needed for quantum state sim¬ 
ulation, it also yields the property that an entangled state cannot be expressed in terms 
of direct product of tensors by using PPS properties and phase ensemble average. In the 
following we use density matrix to illustrate this feature. We assume that a simple type of 
simulation state of N fields can be expressed: 


N' 


N 


Id') = iV'i) 0 |*/’ 2 ) • • • 0 |'0iv) = e' 


iA® 




aO) 


Xj) 


( 6 ) 


2=1 




where N' + N” = 2^ and N” < N, = |h ^2 • • • ^n)- A density 

matrix p can be calculated: 


N' 


N 


N' 


N 


p= 1^) (T| =e 




Ci \xi) + Y^ ^ l^i) 1 






*„-iAO) 




2=1 


J = 1 


2=1 




(7) 
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which is simplified into 


N 


P = '^\Cn\^\Xn) {Xn\+ ^ {C*,Ci\Xi) {x^i\ + C*Ci> \x^') {Xi\) 


( 8 ) 


n=l 


iz^i'=\ 


N 

+ {c}Cie'’''-‘'\xi){x.\+qCfe-<>-"'\x.){x,\) 


n' n” 


i=l j=l 


E E \^j) {^i\ + \x^) {x 


where — X^^'\ By applying phase ensemble averaging 23|, mean reduced density 

matrix is defined p = ^ J2k=i P- balance property of BPS, then we obtain 


2^ n' 

P='^\Cn\^\Xn) {Xn\+ {C*, Ci \Xi) {x^'\ + C* Ci' Ix^') {Xi\) (9) 

n=l 

Eq. on]) shows that all non-diagonal terms including \xj) disappear and the reduced density 
matrix p might not be expressed in terms of a direct product of the states |x„), similar to 
the case of quantum entanglement states. 

In addition to the fact that a quantum entanglement cannot be expressed in terms of direct 
tensor product, quantum entanglement also make a correlation measurement different. The 
correlation analysis on the simulation states is necessary because the nonlocal correlation 
with Bell’s inequality and equality criterion is the most fundamental property of quantum 
entanglement. In order to perform the correlation analysis, a correlation measurement P on 
l"^) is given 


F(0) = q|P(0)|q 


a* P* 



+ af5*e-^^ 


( 10 ) 


For convenience, coefficients a, 13 are set to be l/\/2, yielding P{9) = cos(6*). Further we 
generalize P to the case of N fields: 


P{9i, ...9 n) = P (0i) 0 P {92) 0 ... P (0iv) (11) 

Then we obtain the correlation analysis of the simulation states using P and the density 
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matrix p\ 


N 




= Tr 


pf>(9i,...9jv) 


( 12 ) 


k=l 


y~! |C'n 


N' 


{Xr. 


n=l 


p\xn)+ Y1 (c*>Ci{x,>\P\x,) + c:apxi\P\x^>) 


Eq. ffT^ shows that only non-diagonal terms J2^i'=i I ^ 1^*) + (^*1 ^ l^i') 

remain. 

Key to an effective simulation of quantum entanglement is that the physical resources 
for the simulation does not increase exponentially with number of particle. In the follow¬ 
ing we discuss analysis of computation complexity. A simple unitary transformation, NOT 
gate, is used as an example to show computation complexity. Starting with a single held 
\Pn) = (an | 0 ) -|- f3n | 1 )), applying a unitary transformation switching 17 : | 0 ) ■(-)■ | 1 ) to 
decomposes PPS into each phase unit: U \'ijjn) U («„ |0) -|- /?„ |1)) |fc = 1... 77 = 

(an |1) + f^n |0)). For each phase unit, its computation is the same as that in quantum 
computation, therefore computation for N phase units equals N times of quantum compu¬ 
tation of each phase unit. We can extand unitary transformations to simulation states with 
N helds: 


u 

I'h) ^ 

1^') 





1 

1 

■ N 


- 

1^) = 


E 


A . 

. 4 • • • "In) 


il=0 

ijv=0 

-i=i 


- 


1 

1 

■ N 


- 



E 

rOY 

ii. 

• 4... 


il=0 

ijv=0 

-i=i 


_ 


and coefficients and are related by an unitary transformation: 


(13) 


C 


UY 


*1- 


■In 


E Tjik(~<U) 




(14) 


U) 

Because a PPS contains N phase units and N time slots, therefore the required com- 

ri 

putation is N times that of quantum computation, but 2^ times is unnecessary [3j. 

Two-particles Bell states: Consider the case that the modes |1) in the states \'ipa) and 
I'i/'fe) similar to Eq. ([T]) are exchanged by a mode exchanger constituted by mode splitters 
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and combiners 


2ll . |26| . The exchange yields the following states: 


V2 ^ 


(15) 


V2 


where the relative phase seqnences (RPSs) 7 *-“^ = — 7 ^^) = — X^°‘\ and 7 *-“^ + 7 *^^^ = 0 

The simnlation state |T) is obtained: 


|T) = 


A)® 


A) = 


|0)|0) + |l)|l)+e‘’'‘’|l)|0) + e‘’'‘’|0)|l) (16) 


Appearently the rednced density matrix p cannot be direct prodnct decomposited dne to 
only non-diagonal term | 00 ) ( 11 | -|- | 11 ) ( 00 | remaining. 

Then we obtain the resnlts of the helds in the correlation measnrement P{6a, k) = cos(0^-7 
7 ^,“^) and P{9b,k) = cos{9b -|- where 7 fc“\ 7 fc^^ are the k-th nnits of the RPSs 7 ^“^ and 
respectively. Then the correlation fnnction is 

1 ^ 

E{9a, Ob) = ( 1 '^) 

k=l 


where C = 1/2 is the normalization coefficient. The states in Eq. flT^ are considered to 
be a classical held simnlation of the Bell state jd'/A By snbstitnting the above correlation 
fnnctions into Bell ineqnality (CHSH ineqnality) 27| : 


B\ = |-B(9„,9i) - + £(91.9;) + £(91.9i)| = 2^2 > 2 


(18) 


where 9a, 9'^, 9^ and 9'^^ are vr/d, —vr/d, 0 and 7r/2, respectively, when Bell’s ineqnality is max¬ 
imally violated. 

Bell state |T+) differs from |d/~) by tt phase. Similarly, simnlation of the Bell state |T~) 
is expressed as \Aa) = (^|0) + |l)j /\/2, \Ah) = (^|0) -|-|l)j /\/2. 

By performing the transformation '■ |0) 77 |1) on of th® state 

|d/^), we obtain the simnlation of the Bell state expressed as \'ip'a) = 

^|0)-t-6*'*'''“'|l)j/\/2, j-i/fe) = ^|1)-t-|0)j/\/2, and of |$“) expressed as 

\'^'a) = (^|0) + P^^'"'’ |l)j /\/2, \i/j'b) = ^|1) -1-|0)j /\/2. Then their cor¬ 

relation fnnctions E^,- {9a, 9b) = — cos {9a + 9b), E^± {9a, 9b) = ± cos — 9b) are obtained. 
To snbstitnte the correlation fnnctions into Eq. flT^ . we also obtain the maximal violation 
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of Bell’s inequality. The violation of Bell’s criterion demonstrates the nonlocal correlation 
of the two classical helds in our simulation, which results from shared randomness of the 
BBSs. 


GHZ states: The nonlocality of the multipartite entangled GHZ states can in principle be 
manifest in a new criterion and need not be statistical as the violation of Bell inequality 28 |. 
Breparing three states I'lpa), iV’fe) ^md \'ipc) similar to Eq. ([1]), and by cyclically exchanging 
the modes | 1 ) of the states, we obtain the states as following 


^:) = ^(io)+£,‘v>|i)) (19) 

= ^(|0)+e‘^'*’|l)) 

d) = ^(|0>+e‘’“’|l)) 

where the BBSs 7 *^“^ = 7 ^^^ = A^'’^ — A*^^\ 7 *^'’^ = A*-®^ — and 7 ^“^ + 7 ^^^ + 7 ^'’^ = 

0. We obtain the measurement results P{6aik) = cos{6a + 7 ^“^), P( 6 *fe,/c) = cos( 6 'b + 
lk^)^P{^c,k) = cos( 0 c + 7 f^) for the states \'ip'a) ,\^'b) and in the correlation mea¬ 
surement, respectively, and the correlation function 

1 ^ 

E{9a, Ob, 0c) = j;^Yl ( 20 ) 

k=l 

where C = 1/4 is the normalized coefficient. If 9a + 9b + Oc = 0, E{9a,9b,9c) = 1. If 
9a + 9b + 9c = vr, E{9a, 9b, 9c) = —1. By using GHZ State, the family of simple proofs of 
Bell’s theorem without inequalities can be obtained [26], which is different from the criterion 
of GHSH inequality. The sign of the correlation function can be also treated as the criterion, 
such as the negative correlation for nonlocal and the positive correlation for local when 
9a = 7i/3,9b = 7 i/3,9c = 7i/3. We obtain that the simulation state in Eq. (IT^ shows the 
negative correlation. The results are similar to the quantum case of GHZ states. 

Further, the simulation of GHZ state could be generalized to the case of N particles. By 

preparing N states similar to Eq. ([1]) and cyclically exchanging the modes |1) of the states, 

the RBSs satisfy 7 ^^^ 7 *-^^ = 0. We obtain the correlation function 

1 ^ 

E{9i, ■■■0 n) = j;^Y1 ^(^1’ = cos (01 + ■ ■ ■ + 9n) (21) 


where P{9i, k) = cos(0j 4 - 7 ^) is the result of the classical field with i-th RBSs at the fc-th 
sequence units in the correlation measurement, and C = is the normalized coefficient. 









Using the same notion, we can obtain simulation results of other quantum entanglement 
states. It should be pointed out that the phase randomness provided by PPSs is different 
from the case of quantum mixed states. Quantum mixed states result from decoherence and 
all coherent superposition items disappear. In contract to the decoherence, some coherent 
superposition items remain in the simulation state due to the constraints of the RPSs, such 
as + 7 ^^^ = 0 , 7 *^“^ + 7 ^^^ + 7 *^^^^ = 0 for the simulation of Bell states and GHZ state, 
respectively. These remaining items make it possible to simulate quantum entangled pure 
states. 

As shown in the above examples, we utilized the properties of PPSs to label classical helds 
that are even overlapped in the same space and time. In simulation of entangled states, the 
resources required are the PPSs instead of classical held modes. It means that the amount 
of PPSs grows linearly with the number of quantum particles. According to the m-sequence 
theory, the number of PPSs in the set equals to the length of sequences, which means that 
the time resource (the length of sequence) required also grows linearly with the number of 
the particles. 

In conclusion, a novel simulation method for quantum entanglement is presented, with its 
mathematical expressions and physical meanings identical to those in quantum mechanics. 
In the framework of quantum mechanics, the overall phase of a wavefunction can be ignored, 
as it has no contribution to the probability distribution. However, quantum entanglement 
must be related to two or more spatially separable quantum particles. By introducing a phase 
factor to superposed states with PPS properties, we conclude that quantum entanglement 
can be efficiently simulated by using a classical held modulated with PPSs. The research 
on this simulation not only provides useful insights into fundamental features of quantum 
entanglement, but also yields new insights into quantum computation. 
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Fig.l Construction pathway of simulation states is displayed, which are generated by uni¬ 
tary transformed from initial states-the mode superposition of classical helds with 

PPS A«. 
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